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Abstract 

We analyze a two dimensional SU(3) gauge model of Wil- 
son lines as a dimensionally reduced model of high temperature 
QCD^. In contrast to perturbative dimensional reduction it has 
an explicit global Z3 symmetry in the action. The phase diagram 
of the model is studied in the space of two free parameters used to 
describe the self interaction of the Wilson lines. In addition to the 
confinement-deconfinement transition, the model also exhibits a 
new Zs-breaking phase. These findings are obtained by numer- 
ical simulations, and supported by a perturbative calculation to 
one loop. A screening mass from Polyakov loop correlations is 
calculated numerically. It matches the known QCD3 mass in a 
domain of parameters belonging to the normal deconfined phase. 
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1 Introduction 



It would be of considerable interest to understand the long distance properties 
in the quark gluon plasma quantitatively. Direct simulations of QCD are still 
too costly to be continued to the continuum limit, at least as far as these 
properties are concerned. Dimensional reduction has been shown to be a 
reliable approximation, if the temperature is high enough. In this approach, 
the non static modes of the fields are integrated out perturbatively, in practice 
usually to one or two loop order. The resulting effective theory of the static 
mode is a bosonic field theory in one dimension less PP~[S]- This becomes 
a very powerful technique, when the long distance properties are calculated 
non perturbatively on the lattice in the much simpler reduced theory |B]-|16j. 
In general one finds that if the temperature is larger than about twice Tc, 
the agreement with the full theory in pure gauge theories, where it can be 
checked, is good. A review of the situation can be found in [T^ . 

In a series of articles ITHl IT^ we have investigated the dimensional 
reduction of SU(3) gauge theory in three dimensions to a two dimensional 
Higgs-gauge-model. In this case one can also easily simulate the full theory 
and compare the results in detail. We found that using the one loop approx- 
imation in the perturbative integration of the non static modes, the reduced 
model describes very well the correlations between Polyakov loops as well 
as spacelike Wilson loops down to about 1.5 times the critical temperature. 
It should be emphazised that no free parameters are involved, because the 
lowest order terms of the self-interaction (potential) of the Higgs field are 
determined by the perturbative calculation. 

However, this reduced model does not have the Z3 symmetry of the full 
model. It therefore does not have the confinement-deconfinement phase tran- 
sition. As we pointed out in [20^ one may construct a model, which has the 
perturbatively reduced model as the high temperature limit, but which does 
not break the Z3 symmetry. Similar models have been discussed in |22] • 

In this paper we investigate a model where the basic degrees of freedom 
are the spatial components of the gauge field and the Wilson lines of the three 
dimensional pure gauge theory. The traces of the latter are the Polyakov 
loops, order parameters for the preserved Zs-symmetry. The Polyakov loop 
potential is left undetermined and we parametrize it by the first two terms 
of a small loop expansion. We study the phase diagram of this model in the 
corresponding parameter plane and show that in addition to the expected 
confined and deconfined phases, another Zs-breaking phase, distinct from 
the normal deconfined phase, appears. 

The Polyakov loops correlation length is also calculated at high tempera- 
ture along the line corresponding to a purely quadratic potential. It matches 
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that of the original 3D gauge theory at a point which hes inside the normal 
deconfined phase, not the new one. In this region of parameters, it is thus 
consistent to view the model as a dimensionally reduced model for QCD3, 
although an effective gain in computing cost requires an a priori knowledge 
of the Polyakov loops potential, not calculated in this paper. 

These calculations are performed via Monte Carlo simulations. First ac- 
counts of the numerical results have been presented in |231I21I- Furthermore 
we show analytically that the one loop approximation gives a qualitative 
description of the phase diagram. Although our numerical simulations are 
restricted to the two dimensional reduced model, the one loop approxima- 
tion can also be directly applied to the three dimensional case, which is the 
reduced model of the physical four dimensional QCD. We see in this approx- 
imation the same qualitative behaviour of the reduced model, as in the two 
dimensional case, including the existence of the new phase. 

In section 2 we define the reduced model, in section 3 we describe the 
numerical simulations, and in section 4 we give the corresponding results. In 
section 5 we derive the effective action of the reduced model in the one loop 
approximation, and compare with the results of the simulations. We also 
present the analytic result in the corresponding three dimensional reduced 
model. Finally section 6 contains a summary of our results. 

2 The Model and its Relationship to Finite 
Temperature QCD^ 

In this section, we first write down the action of the 2D-model under study 
and then discuss under which conditions it can be viewed as a dimensionally 
reduced action for the Polyakov loops of pure QCD in 3 dimensions. 

The model is defined by its partition function on a square lattice with 
sites X, spanned by unit vectors i, i=l,2. The spacing a will be usually set to 
1, unless necessary when dimensionful quantities are defined. All dynamical 
variables belongs to the SU(3) group, gauge fields t/(x, i) sitting on the links, 
and matter fields ^(x), localized on the sites. The variables ^(x) will be 
later identified with the Wilson lines of QCD3, whose normalized traces are 
the Polyakov loops. The partition function is written 
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VUVV exp {-Su - Suy - ^y), (1) 
Su = /32 5^('l-i3ftTr[f/(x;l)f/(x+T;2)f/t(x + 2;l)f/t(x;2)]Y(2) 

X i=l,2 ^ ^ 



(3) 



The measure VU W represents the products of the Haar measures asso- 
ciated with the U and V variables, Su is the 2D Wilson action and Sjjy a 
gauge invariant kinetic term for the matter fields, defining lattice couplings 
P2 and t. The self-interaction Sy is assumed to be a real and locally gauge 
invariant function of V, to be specified later. 

Actions of the above type have already been considered in El I22j 
as representative of what an effective action for the Wilson lines of QCD^ 
might look like. A first account of some of our results has been presented in 



The connection of the partition function to QCD^ may be thought of 
as follows. Adding a third dimension (a time direction 0), of length aLo, to 
the previous lattice, QCD^ at temperature T = l/(aLo) is defined by the 
3-dimensional Wilson action on this new lattice and a coupling (3^: 

Lo 2 

Sqcds = /^s E E E 

xo=l X ij,<u=0 

1 - ^3? Tr U{x; fi) U{x + Jl; z/) U{x + u; fi)^ U{x; u)^^ . (4) 
The continuum limit is the limit a — when 

9l = (5) 
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which sets the energy scale, is kept fixed. So we define a dimensionless 
temperature by 

T (3, 



(6) 



and fixing the temperature amounts to letting (3^ and Lq go to 00 in a fixed 
ratio. Of course at the same time, the spatial lattice size a L is supposed to 
be large compared to any correlation length. Dimensional reduction Ej 
provides a way to extract properties of QCD in the large T / gl regime. The 
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application of the general procedure OIHEl to QCD^ is described in details 
in m- 

Here we only recall the main steps of the process. In the Wilson action 
(jH), the dynamical field variables U{x;fi) G SU{3) are related to the gauge 
fields by 

U{x;fi)=exp[iA^{x)]. (7) 
To make the reduction, we choose a static gauge 

^o(a:o,x) = Ao(x), (8) 

which means 

f/(xo,x;0) = t/(x;0) (9) 

for all xq. The variables V{x.), which we call Wilson lines, are defined in this 
gauge by 

Vi^) = f/^"(x;0) = exp[zLo Ao(x)], (10) 

and are static operators. Their normalized traces are the gauge invariant 
Polyakov loops, denoted f (x): 

^;(x) = ^Trl^(x). (11) 

The spacelike fields Ai{xo, x), i = 1, 2 , can be split into static and non-static 
components in this gauge, 

^•(x) = f J]A,(xo,x), (12) 

xo 

Ar{xo,x) = A(xo,x)-A^(x). (13) 

The last step of standard dimensional reduction consists in integrating per- 
turbatively over the A^'^ fields. For this purpose, the lattice action (jH) is 
expanded to an appropriate order in powers of Aq and Ai, and full gauge 
fixing is achieved |25j \2M setting the condition 



(14) 



( Aj(xo,x) - Ai(xo,x-i) ) = 0. 

Together with Eqs. ()12ll3p . it implies 

E (^•(x)-AI(x-T)) =0, (15) 
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which is a lattice Landau gauge for the 2D gauge fields (x). By integrating 
perturbatively over the gauge fields A^'^{xo,x.) up to one loop order, one 
gets a two dimensional effective gauge theory of the static modes, where 
A|(x), i = 1,2, are the gauge fields and Ao(x) acts like a Higgs field in the 
adjoint representation. There are no infrared divergences in this perturbative 
calculation, and one obtains at a given order in (yf|/T only a finite number of 
interaction terms. The effective two dimensional model has, however, severe 
infrared divergences in perturbation theory and thus has to be treated non 
perturbatively. The lattice version of the reduced action is 



S2L — Su + Su^Ao + Sao: (16) 
Su = P3L0 5^ (1 - Tr t/(x; l)t/(x + T; 2)t/t(x + 2; l)t/t(x; 2)^ , 
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^EETr(A(f/)Ao(x)) , (17) 

X j=l ^ ^ 

A(^7)Ao(x) = f/(x;2)Ao(x+?)f/t(x;^) - Ao(x), 

Sao = Y.^2 Tr A2(x) + (Tr A2(x))' . (18) 

X 

In the above formulae, Su is the pure gauge term in 2D, Su^Ao the gauge 
invariant kinetic term for the scalar Higgs field and Saq the scalar potential, 
whose self couplings /c2 and /C4 result from the one loop integration over the 
non-static components of the gauge fields ^Hj- All terms have a global R-r 
symmetry Ao(x) — y4o(x), while the symmetry of the original 2+1 
dimensional SU(3) is broken by the perturbative reduction procedure. 

As shown in JH] ■, the 2D model defined by the lattice action ()16|) success- 
fully accounts for non-perturbative properties of QCD^, at T > 1.5 Tc and 
large distances, a conclusion reached from the behaviour of Polyakov loops 
correlators. 

We now examine under which conditions the model of Eq.(m may also 
be related to this perturbatively reduced model. First the pure gauge piece 
Su of Eq. (j21) coincides with that of S2L if 

P2=LoP3. (19) 

Next consider the content of Su,v, Eq.Q, for V given by (fTO|l and Aq small, 
the situation in which dimensional reduction is known to be justified. To 
second order in Aq, Eq.® generates the lattice invariant kinetic term 



X j=l,2 



^^(x) + Ag(x +T) - 2 t/(x; i)v4o(x +T)f/^(x; i)Ao(x) 



(20) 
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which coincides with that given by Eq. (fT7j) provided one chooses 



t = Ps/Lo. (21) 

In other words, the couphng t of the model is proportional to the dimension- 
less temperature defined by Eq. (jH)). 

We finally come to the last term in ((H), unspecified yet. Of course one 
would like to be able to compute it from QCD^, as it is the case for the self 
interaction 5*^0, whose lowest order terms (fTKjl were computed perturbatively 
in the high T and large distances situation 16^. In the absence of a known 
scheme to do so, we choose the simple form 

Sv = J2 (^2 \v{^)\' + As ^v%^)] , (22) 

which fulfills the following requirements. It is real and locally gauge invariant. 
It is also invariant under the global transformation 

V(x) ^ ZnV{^), (23) 

Zn = exp(2i7rn/3); n = 0,1,2 mod{3). 

This in fact true for the full action of our model, and distinguishes it from 
perturbative dimensional reduction. Being an SU(3) matrix, l^(x) possesses 
two (invariant) degrees of freedom only, for example two of its three eigen- 
values 



K,(x) = exp 



(24) 



n = 1' (25) 



a=l 



or equivalently the real and imaginary parts of its normalized trace w(x), an 
order parameter for the Z3 symmetry. A self-interaction like ()22p opens the 
possibility of a transition from a symmetric phase, where (f (x)) vanishes, 
to a broken phase where it has a finite modulus. From this point of view, 
we consider fl22|l as the first terms of an expansion of Sy near the phase 
transition. The attempt made in [21]] was to fix A2 and A3 in such a way 
that, for Aq small, Sy matches the two lowest order terms of the expansion 
(dl) of 

The use of free parameters in Sy makes the model less predictive than 
perturbative reduction at high temperature. In turn, a parametrization like 
()22|1 allows to investigate which non-perturbative effects might be associated 
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with lower temperatures. This is why, in what follows, we focuss on the 
structure of the phase diagram, which we study numerically in the 1X2,^3} 
plane at A3 = 0, the effect of a third order term in v being discussed at a 
perturbative level only (section 5). 

3 Numerical simulations 

The model (0) was studied using conventional lattice QCD Monte-Carlo tech- 
niques. Lattices of sizes ranging from = 16^ to 72^ were studied and the 
parameter Lq (the inverse temperature of the 2 + 1 model in lattice units) 
was always set to 4. The multi hit metropolis algorithm (with 8 hits) was 
used for updating both the gauge fields f/(x, ?) and the Wilson lines V^(x). 
The dynamics of the two sectors were quite different: While the integrated 
autocorrelation time for the plaquette operator was negligible (of the order 
of 10 sweeps), for the Wilson lines it could rise up to 25000 sweeps near the 
phase transition on a 32^ lattice. One sweep consisted of one update of all 
the U fields followed by one update of all the V fields. Typically, runs with 
2-4 10^ sweeps were performed for the system on a 32^ lattice around the 
phase transitions. 

The quantities measured concerned the Polyakov loops f (x) defined in 
Eq. (ITT|) : their lattice averages 

X 

their two-point on axis correlators 

= ;^($^$^^^(a;i,X2)t;*(xi + r,a;2)) (27) 

Xl X2 

and projected correlators : 

Xl XI x'2 

The traces L were measured and stored on a per configuration basis, so 
that their distribution P(£) and the susceptibility 

x = m^)-m? (29) 

could be obtained. 



8 



The correlators were fitted with the appropriate formulas (see section 4.2) 
using correlated fit procedures in order to extract screening masses fI7\ . 

All the errors were calculated using a blocked bootstrap method. The 
data were divided into blocks, the size of which was adjusted so as to make 
them independent from each other, but at least ten blocks were used. Then 
a bootstrap sample was drawn from those blocks and the measurements per- 
formed on this sample. The procedure was repeated 200 times and the stan- 
dard deviation of the resulting distribution of measured values taken as an 
estimate of the error on their average. 

4 Results from Numerical Simulations 

The numerical simulations gave us evidence for the phase diagram in Fig. ^ 
as we explain below in 4.1. Further in 4.2, we discuss the behaviour of the 
correlation length measured in the various phases. 

4.1 Phase structure 

We started with simulations of the so-called 'naive' model, i.e. with no self- 
interaction of the Polyakov loops: Sy = in Eq.((T)), that is A2 = A3 = with 
the parametrization (j221)- There the model may also be viewed as standard 
2-1-1 QCD on a lattice with one time slice only. Exploring a range of /?3 values 
we observed a clear signal for the phase transition expected between a low 
temperature (low Ps) confined Z[3) symmetric phase (region la in Fig.Q) and 
a broken symmetry high temperature phase (region II). The corresponding 
peak in the susceptibility occurs at very close to 20 as shown in Fig. |2l 
where the average of |£| is also shown. 

Still keeping A3 = 0, we then performed series of scans in the /53-A2 
plane on 16^ lattices, looking for peaks in the susceptibility as signals 
for phase transitions and characterizing different phases by the corresponding 
distributions of C. Although it also includes informations from larger lattices, 
the tentative phase diagram of Fig. ^ results from this exploration. 

More precise scans in A2 along the line (3^ = 42 (T ~ 2Tc) give for the 
susceptibility the results plotted in Fig. El A strong transition signal shows 
up just above A2 = 0.3. There, the time history of |£p and the corresponding 
histogram (Fig. clearly favour a first order phase transition. The height of 
the peak is furthermore compatible with a size dependence in the range 
L = 24-48, which is the finite size scaling expected for a first order transition. 
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Figure 1: A tentative phase diagram. Open points denote the location of 
peaks in susceptibihty observed on 16^ lattices. Black points denote the 
position of phase transitions established on 32^ lattices with large statistics. 
The lines are hand drawn sketches of where transitions take place. 
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Figure 2: Polyakov loop susceptibility (circles) and average modulus (dia- 
monds) versus Ps for A2 = A3 = on 32^ lattices. The continuous lines result 
from Ferrenberg-Swendsen reweighting. 
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Figure 3: Polyakov loop susceptibility versus A2 for = 42 and A3 = 
on lattices of various sizes. The continuous lines result from Ferrenberg- 
Swendsen reweightings for the three smallest lattices. 





20000 40000 60000 80000 0,1 0,2 0,3 0,4 0,5 



Figure 4: Time history and histogram for the quantity \jC\'^, on a 32^ lattice 
for — 42, X2 = 0.305 and A3 = 0. One clearly sees jumps between two dis- 
tinct values, and the corresponding double peak structure of the histogram. 
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The evidence, from similar considerations, for a second phase transition 
around A2 = 0.4 is considerably weaker. A peak is hardly visible on the 
same lattices and we had to go to 72^ lattices (with correspondingly poorer 
statistics) to ensure that the bump detected does grow with the size. As 
we will immediately see however, the combination of the above observation 
with properties of the Polyakov loop distributions confirms the existence of 
a second transition. Furthermore, the theoretical analysis developed in the 
next section supports the same conclusion. 

In order to identify the nature of the various phases expected from the 
diagram of Fig. we looked at the corresponding distributions of Polyakov 
loops. The phases labelled la and II are easily identified with standard 
symmetric and broken Z[3) phases respectively. In Fig.|21 these distributions 
are shown in the complex plane of C, and they exhibit the expected patterns: 
The values of C cluster around in phase la, and have a finite modulus and 
an argument of the form 2i7rn/3 in phase II. 

As we increase A2 from at /^a = 42 and move from region II to region III, 
another symmetry breaking pattern appears. symmetry breaking is again 
manifested by clusters around non-zero values of C, but now arguments of 
the form 2iTT {n + l/2)/3 are favoured. This is illustrated on Fig. IHl (left). 
In this example, the lattice is not large enough to prevent the system from 
tunneling between different phases. Also the run is too short to allow for full 
equilibration. But accumulations of values of C close to the boundary with 
arguments around vr and — 7r/3 are clearly observed. Increasing further A2, 
we move into a region lb where Zs-symmetry is recovered, as shown in Fig. El 
(right). We conjecture that this phase is the same as the low temperature 
confined phase (la). 

4.2 Screening masses 

The inverse of the correlation length between Polyakov loops in the present 
model corresponds to a screening mass in 2+1 QCD. This screening mass 
was obtained from the two point Polyakov loop correlators ^T7\i by fitting 
m5 in 



The constant c accounts for the fact that we measured unconnected correla- 
tors. Fig. 13 provides a comparison of ms, as obtained in the present model 
for A2 = A3 = (circles), with previous results (black squares) obtained in 
the naive perturbative reduction scheme {^^, no Higgs potential). The data 
are plotted versus the quantity T/T^ = Ps/Psc, where is the critical value 




(30) 
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Figure 5: Typical distributions of £ at A2 = A3 = for phases la (/^s = 18.0, 
left) and II {P3 = 24, right). 




Figure 6: Typical distributions of C cd, (3^ — 42 and A3 = for phases III 
(A2 = 0.32, left) and lb (A2 = 0.48, right). 
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Figure 7: Screening masses in lattice units for the model (P) with 5*^ = 
(circles). At large T, they approach those obtained in naive dimensional 
reduction of QCD^ (black squares, 2D adjoint Higgs model with no potential 

m)- 

of /^s for Lq = 4, respectively equal to 20 ( this model, see Fig. H]) and to 
14.7 in 2+1 QCD j2H]- The two sets of points tend to join each other at 
high temperature. It is the expected result since if V becomes equivalent to 
an element of Z3, it can be rotated to the form ()1()|1 with Aq small, and the 
action in (^J) can thus be expanded in powers of Aq. In this situation, if there 
is no loop potential Sy to generate a local Higgs potential, one ends up with 
the naively reduced model. Near the transition, the two models have very 
different behaviours. In the perturbative reduction there is no symmetry 
in the action, correspondingly no phase transition, and the mass does not de- 
crease near the (5^c of the full theory. In the present model, there is a phase 
transition and the mass decrease from above the transition is consistent with 
a second order transition. This in accord with the behaviour of QCD^. 

We also measured the screening masses as a function of A2 from our data 
at /^s = 42 and A3 = 0. In this case, we used the projected (0-momentum) 
correlators (|^Hj) and extract the masses by fitting 777,5 in 

Pp,^(r) ^ a(e-™s- + e^"^^^^-")) + c. (31) 
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Figure 8: The screening masses ms in lattice units, fitted with Eq. at 
(3^ = 42 and A3 = 0, as a function of A2 for various lattice sizes. 

Our results are summarized in Fig. |Hl for three lattice sizes. Important size 
effects are observed in the region of the (/// —>■ lb) transition just above 
A2 = 0.4. It is not so near the transition at A2 = 0.3, which is consistent 
with evidences mentioned earlier which favour a first order transition. It 
is important to notice that the screening mass at (3^=4:2 (corresponding to 
T/Tc ~ 2) coincides with the screening mass of QCD^ at A2 = 0.18, that is 
well inside phase II. The situation is different from that in the perturbatively 
reduced model, where we found the reduction point to be in a metastable Z2 
symmetric region of parameter space, beyond the transition to the Higgs 
phase UHlini. 

5 Perturbative Approach to the Polyakov Loops 
Potential 

In order to get further insight into the phase structure of the model, we 
have performed a perturbative calculation of the effective potential of the 
Polyakov loops. We assume that in a phase where Z3 is spontaneously bro- 
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ken, i.e. where (f (x)) is not zero, the Wilson hne ^(x) is, up to a local 
gauge transformation r(x), a constant matrix which we denote V. So any 
local fluctuation of V^(x) which cannot be gauged away is ignored, and the 
Polyakov loop f (x) is a constant number v on the lattice: 

v(x) = r(x)vrt(x), (32) 

i;(x) = ^TtV{^) = \ttV = v. (33) 

Due to gauge invariance of both the action and the measure, r(x) can then 
be absorbed in a redefinition of the group variables U. According to Eqs. 
((131221), the Polyakov loops potential Wp{v) defined by 



exp[-L^Wp{v)] = j D[f/]exp(-5), 



IS 



Wp{v) = \2\v\^ + X3'Slv'^ + W{v), (34) 

where W{y) is the contribution resulting from the coupling of V to the gauge 
fields. In order to compute it to one loop order, we first expand Sjj + Sjjy 
to second order in the SU(3) algebra using 

f/(x;i) = exp[i5( Ai(x)], (35) 

8 

A,(x) = 5^r'^4«(x), 

a=l 

Tr[T'^T^] = ]^6ab. 

Note that with respect to previous definitions, Aj(x) has been rescaled by 
the gauge coupling g, where = Q/P2- Including the quadratic part of the 
jacobian from the U A change of variables, we get 
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J2 ^?(x) - Ai(x + 2) - /l^(x) + A^(x + 1) 

. a=l \ 

+ E ( 4^ E ^"(^)^' W (^^at - M., - M,^ + ^ (A^(x))^ 

i=l \ a,b=l ^ ^ a=l 



(36) 



The 8x8 constant matrix M is defined by 

M^b = 2 Tr [r'^1/r*0], (37) 
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and its properties are derived in the appendix. Since integration over the 
Ai's requires gauge fixing, we then add up to S\ the i^-gauge fixing term 

^« = 7 E E (^1 W - ^1 - 1) + ^2(x) - A^(x - 2)] (38) 

Note that the hmit ^ ^ reproduces the Landau gauge (fT^ in which 
the perturbative reduction was performed. There is no contribution from 
the ghost terms to the one loop potential. For the rest of the calculation, 
whose details are given in the appendix, we go to Fourier space, using lattice 
momentum variables 

'^TTTi' ^ ^ ^ ^ 

ki = — ^, ki = 2 sin ki/2, k^ = kj + kl (39) 
Ij 

The quadratic form in Ai, diagonal in fc-space, is a matrix Q{k) in the 
coordinate (3) colour space so that 

^(^) = ^ Y.^ogdetQ{k) = ^ ^TrlogQ(A;). (40) 

kik2 kik2 

Here, Tr denotes the trace over both the coordinate and color indices. We 
obtain 

Qt!; = K.,,5ab-l^Ms,at6i,, (41) 

7' = Ig't, (42) 
K,, = (fc^ + 7^ + ^)5,,-(l-l)^,^,exp(z^^i^), (43) 

Ms = ^(M + M^). (44) 

The diagonalization of K is trivial. That of Ms is made in the appendix, 
where in terms of the eigenvalues (j21I) of V, its eigenvalues are shown to 
be 

1, cos(ei - 62), cos(e2 - ^^3), cos(^^3 - ^^i), (45) 

each of them being doubly degenerate. Up to an additive constant in v, the 
one loop contribution to the Polyakov loops effective potential then reads 

^(-) = ^E E logR(fc^ + ^ + 2f sin^^^) X 

fcifco a<0=l \ / 



where 



72 I c ^9'^ I n 2 • 2 — 6f3- 
k + t{ h 2 7 sm 



(46) 
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Because the eigenvalues exp{i6a) of V are not independent, it is more con- 
venient to express W explicitly as a function of the Polyakov loop v, which 
is furthermore an order parameter for the Z3 symmetry. This again is per- 
formed in the appendix. Specializing to the Landau gauge by taking the 
limit ^ — > and dropping terms not depending on v, we obtain 



Using this expression, we studied the shape in v of Wp{v) for = 42, Lq = 4, 
A3 = looking for its minimum as a function of A2. 

For A2 = this minimum is reached at |f | = 1, i.e. v = exp(22n7r/3), the 
standard Za-breaking phase II. Of course the situation remains unchanged 
for any negative A2, which enhances the Boltzmann weight. On the contrary, 
large 1^1 values are suppressed for A2 positive, and the minimum is at f = 
for A2 large enough, this is phase lb. There is however an interval in A2 
between ~ 0.25 and 0.36 where the minimum occurs at the non trivial points 
f = 1/3 exp[2i (n + 1/2) 7r/3], corresponding to the phase III discovered in 
our numerical simulation inside a similar interval (A2 between ~ 0.3 and 0.4, 
see Fig. 0}. The existence of this intermediate phase is due to the cubic term 
inside the log of W{v) in Eq. which is minimum for arg(f^) = vr. 

The situation just described is illustrated in Fig. IHl On the left, we 
show Wp{v) in the complex f-plane for A2 = 0.32 inside the phase III, to 
be compared with the similar plot of Fig. IHl (left), the Boltzmann weight 
being the largest (resp. smallest) in the blackiest (resp. clearest) regions. 
On the right, Wp{v) is plotted versus v real, for A3 = and three values of 
A2 corresponding to phases II, III, and lb, as seen from the locations of the 
minimum, respectively at v=l, f=-l/3 and v=0. 

The analysis can be extended to include the effect of the term X^^v^. 
From the discussion above, it is clear that a positive (resp. negative) value 
of A3 favours (resp. disfavours) phase III, which is characterized by 
negative. In fact, computing Wp{v) for f =-1/3, 0, 1 in the A2-A3 plane, the 
perturbative phase diagram is easily obtained: At each point the true phase 
corresponds to the lowest of the three Wp values, a (straight) transition line 
occurs whenever the two lowest coincide, and a triple point when all three 
are equal. The result is shown in Fig. E3 
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Figure 9: = 42, A3 = 0. On the left, the Polyakov loop potential Wp{v) 
in the f-plane at A2 = 0.32. On the right, the shape of Wp{v) for v real at 
three values of A2. The positions of the minimum characterize phases II, III, 
I successively from bottom to top. 

The lattice perturbative calculation leading to Eq. ()4fij) is formally sim- 
ilar to that made for SU(2) in Ref. jSHj- There, gaussian integration was 
performed over both the static and non static components of the QCD4 spa- 
tial gauge fields. As a result of this purely perturbative approach, the Z3 
symmetry could not be spontaneously broken. 

We end up this section with a few remarks. 

Remark 1 

The perturbative calculation of ly(l^) can be easily extended to any di- 
mension D. The only change concerns the kinetic part Kij PHj) whose indices 
now run from 1 to D. In the ^-gauge, this K-matrix has D — 1 eigenvalues 
equal to Ki and one to K2, where Ki and K2 are given by Eq. of the 
appendix for k"^ = YliLi ^1- Then in D-dimensions the perturbative result 
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dHj) reads 



Wp(v) 
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2 
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Using it in dimension 3, we find that the phase structure represented in 
Fig. El persists, with similar values of A2 and A3 on the transition lines if the 
parameters 7^ and are kept the same. 



Remark 2 

Heuristic arguments neglecting any entropic consideration help to under- 
stand the phase diagram found, including the new phase III. The part Suy 
of the action is trivially minimized if the commutator [f/, V] vanishes for any 
f/, which implies V E Z^, i.e. in phase II. Conversely, typically for A2 > 
and A3 = 0, Sv is minimized by t'=0 (phase I) and there is a competition. A 
possible compromise is to require V to commute with all the elements of an 
SU(2) subgroup of SU(3), that is to be proportional to a Z2 subgroup of Z^: 
In diagonal form, V =diag[exp(i6'), exp(2^), exp(— 226*)], up to colour index 
permutations. In this subset |f | is minimum for 9 = (2n + l)7r/3, and there 
f = 1/3 exp[i (2n + l)7r/3] is in phase III. 

Remark 3 

The physical region for v is limited by the condition that two (or three) 
eigenvalues of V are equal. The points characterizing phases II and III belong 
to this boundary, and the measure T) V vanishes there. This contributes a 
logarithmic repulsive potential on the boundary. Its effect is expected to 
repell the values of v in phases II and III slightly inside the physical region, 
and to change the exact position of the transitions lines. 



6 Summary 

We have investigated a model for Wilson lines coupled to a SUiK) gauge 
field in two dimensions. This model was inspired by dimensional reduction 
of three dimensional SU (3) gauge theory. It is a Z^ symmetric Higgs model 
with the Wilson line acting as a Higgs field in the adjoint representation of 
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the group. In perturbative dimensional reduction, where the Z3 symmetry 
is exphcitly broken, the self couplings of the corresponding Higgs field are 
determined by the perturbative integration of the non static modes in the 
original theory. Here we have left the selfcouplings free parameters, to be 
eventually determined from the full model. The price paid is a loss in pre- 
dictive power, but the model allowed us to investigate in detail a non trivial 
phase diagram, which includes the confinement-deconfinement transition of 
QCD^. 

In fact, the model exhibits three different phases. The order parameter 
is the trace of the Wilson line, i.e. the Polyakov loop. At low values of j32 we 
find a confined phase, where the Polyakov loop is distributed around zero. 
Above a transition we find for A2 = a deconfined phase, where the Polyakov 
loop is near to an element of Z3, one of the third roots of unity. Keeping 
fixed in the deconfined phase and letting A2 grow we find a first order phase 
transition into a new phase, where the Polyakov loops are situated near the 
boundary of phase space but with a phase half way in between the third roots 
of unity. Here also, symmetry is spontaneously broken, but the Wilson 
line is near to an element of a Z2 subgroup of SU(3). If we increase A2 even 
more we come into a confined phase again. 

We also have measured the screening masses from the Polyakov loop 
correlations. They become small near the deconfinement phase transition. 
If one uses the screening mass to compare with the full model, in order to 
find the values of the coupling constants which correspond to dimensional 
reduction, one finds it to be in the normal deconfined phase. 

Considering this Wilson line model as a possible dimensionally reduced 
model of three dimensional SU(3) gauge theory, and comparing it with usual 
perturbative dimensional reduction one finds that 

- the model has explicit symmetry in the action, in contrast to the 
perturbatively reduced model, where the symmetry is explicitcly broken. 

- it has a confinement-deconfinement transition in contrast to the pertur- 
batively reduced model 

- at this deconfinement transition, which is in the same universality class 
as full QCD3, the screening masses should go to zero. 

- the reduction point in parameter space seems to lie in the physical phase, 
in contrast to the perturbative reduction, where it is in the metastable region 
inside the Higgs phase 

All these properties makes the Wilson line model an interesting extension 
of dimensional reduction, if one wants to approach the phase transition, while 
perturbative reduction can only be used at temperatures higher than about 
two times the critical temperature. 

There are several interesting questions that can be further investigated. 
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In particular one would like to extend the numerical calculation to three 
dimensions, where one would have a dimensionally reduced model of {3 + l)D 
QCD, i.e. the physical theory. Furthermore one should address the question 
if the new phase plays any role at high temperature or high temperature and 
finite density QCD. It is of course also important to restrain the parameters, 
so as to get quantitative predictions without parameters, as it is the case in 
perturbative dimensional reduction. 
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A Calculation of the Perturbative Polyakov 
Loops Potential 

Here we derive the expressions fj45p of the eigenvalues of the M-matrix, and 
the final expression of the potential as a function of the Polyakov loop. The 
following two trace identities in SU(3) will be used extensively 



Tr(T''XT"r) = UttX TiY Tt{XY)\ 

a ^ ' 

^ Tr(T"X) Tr(T"r) = ]^(ty{XY) - i TrX TrF^. 



(49) 
(50) 



Recalling that TrT*^ = 0, one first shows from its definition (jHTjl that the 
matrix Miy^ is real and fulfills 

M{V^)M{V) = M{V)M{V^) = I, (51) 
MP{V) = M{VP). (52) 
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Hence M is orthogonal. In the matrix Q of Eq. (H^, we factorize the matrix 
K of Eq. call Ki its eigenvalues and write 



TrlogQ = TrlogK + X, 
Trlog/sT = log Kl + log Kl, 



X 



V — Tr [K-P 

p=i ^ 



M{V) + M{V^) 



(53) 



Note that exists for any ^ > 0. The Landau gauge is reached in the 
limit ^ ^ The trace involved in X is the product of K^^ + by the 
trace rrip in color space of (M + M^y, 



a=l 



M{V) + M{V^) 



(54) 



In this expression, we expand (M + M^)^ and apply (j52|l : 



rrir, 



4 E E Tr 

q=0 a, 6=1 



jia yq rpb yjq 



Tr 



yfp-Q j^o. yp-Q 
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Repeated use of EDI) to sum over color indices leads to 



g=0 



^9 y1P-i 



Tr 



1 . 



(56) 



Let Va = exp(i 6*0,), a = 1, 2, 3 be the eigenvalues of V; those of V''' are l/Va 
so that 



3 ''T/ \ 2<J-P 



9=0 



\a,l3=l 



1 



(57) 



In Eqs. ()53j) . the resummation over p gives X, and follows. Up to 
additive constants independent of V, X and W are then given by 



X 



W{V) 



i=l 



n 1 



a</3=l 



7^ 



L2 



fciA:2 j=l 



n (a-- 

a</3=l 



2 



(58) 
(59) 
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Eq. proves that the eigenvalues of the symmetric part of M are 3fJ Va/Vp = 
cos{9a — Op), as stated in (jl^ . The next task is to express as a function 
of Tr V , and for this we have to evaluate the product 

n m 

a<(}=l ^ a p / 

for T = 77(2 Ki). Set P = En=o(-l)" and use 14 1^2 ^3 = 1- We have 
Co = 1 and, c.p. denoting the circular permutations of {1 2 3}, 

ci = (v;2 + F^') + c.p. 

C2 = {v,^ + vi) V, (vi + vi) Vi + c.p. 

C3 = iv,' + vi)ivi + vi)iv.^ + v,'). 

Replacing systematically V^ + V^ by Tr V'^ — V^ whenever a, j3, 7 are different, 
we find 

ci = TiV^ TtV - TrV^, 

C2 = TtV^ TiV + TiV^^, (61) 
C3 = TtV^ Tr02 _ i_ 

Now Tr V'"' for n > 1 is a polynomial of degree n in Tr \^ (for general formulae, 
see Appendix A in [HI)- In particular 

TtV^ = {TrVf - 2 TiV^, 
TrV^ = (TrF)^ - 3|TrV|2 + 3. 

Substituting these expressions into Eqs. (j6T|l and using t> = 1/3 Tr V provide 
the desired result: Up to an additive constant, the perturbative part of the 
Polyakov loop potential is 



L2 



2 ' 2 ' ' '/II' 2 ' ' 23 

where the eigenvalues Ki are 



,(62) 



i^i = fc^ + 7^ + ^ ; i^2 = | + 7^ + ^ (63) 
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